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Abstract. Deviations from the Friedmann-Robertson- Walker (FRW) cosmological model 
in the form of density inhomogeneities induce observational effects on the light propa- 
gating through these fluctuations. Using a rigorously parametrized metric whose pseudo- 
Newtonian potential is related to (possibly nonlinear) density inhomogeneities through a 
relativistic Green function, the behavior of radiation propagating through this approxima- 
tion to our universe is investigated, without relying on any spatial averaging or late time, 
short perturbation wavelength Newtonian limit. In certain regimes the energy shift of a 
photon due to density fluctuations, the Sachs- Wolfe effect, is found to deviate significantly 
from the linearized relativity (Poisson- Newton) result, but the corresponding cosmic mi- 
crowave background temperature fluctuations are still below 10 -6 on small scales. On large 
scales this can be treated as an effective transfer function for the density power spectrum, 
altering the scaling of the amplitudes of large versus small scale power. 
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1. Introduction 

A self consistent treatment of observations in a realistically inhomogeneous universe in- 
volves relativistic cosmological calculation of both the influence of the inhomogeneity on 
the metric and the metric on the light. Often both are treated only within the late time, 
short wavelength Newtonian approximation where the density determines a spatial metric 
fluctuation by the Poisson equation and the metric influences the propagation through an 
essentially local impulse. Here we use recent improvements beyond "Newtonian" order to 
investigate relativistically light propagation, especially from the last scattering surface of 
the cosmic microwave background radiation (CMB). 

In §2 we use a relativistically rigorous metric (Futamase 1989) describing a universe 
deviating by gravitational potential perturbations from FRW, accurate to a well defined 
order, to obtain the Einstein field equations. The Green function solution of Jacobs, Linder, 
& Wagoner (1992=JLW1; 1993=JLW2) relates the gravitational potential to the density 
inhomogeneity for arbitrary density contrasts, i.e. without restrictions to the linear regime. 
We identify those regions in parameter space where post- Newtonian effects are appreciable, 
as well as deriving analytic expressions for the derivatives of the potential, useful in light 
propagation calculations. In §3 we concentrate on the Sachs- Wolfe effect, where density 
inhomogeneities generate temperature anisotropics in the CMB, evaluating the deviation of 
the Green function results from the standard Newtonian ones. The results are given both 
numerically and by analytic order of magnitude arguments to reveal where the deviations 
arise. In addition we find an effective correction to the density power spectrum which 
alters the intrinsic scaling of its amplitude from small to large scales. 



1 



While the redshift of light propagating through inhomogeneities is well studied (see, 
e.g., Ma & Bertschinger 1995 for the Sachs- Wolfe effect in linear perturbation theory and 
Seljak 1996 for the Rees-Sciama effect in numerical simulations), it has been considered 
within linearized general relativity and not the post-Newtonian formalism of Futamase 
1989, JIW1, and JLW2. That possibly significant differences may arise can be seen from 
the diffusion equation analogy of JLW2 as well as the relativistic approaches of Kodama & 
Sasaki 1984 and Futamase & Schutz 1983, which agree with the results here that there exist 
corrections to the linearized general relativistic solution - the Poisson- Newton equation be- 
tween the gravitational potential perturbation appearing in the metric, 0, and the density 
fluctuation 5p. That is just the late time, short wavelength "Newtonian" approximation. 

2. Gravitational potential 

Under the assumptions that the gravitational potential fluctuations are small, para- 
metrized by e 2 1, and their peculiar accelerations V0 ~ e 2 /n are small (thus ensuring 
that peculiar motions remain much less than the speed of light) for characteristic inho- 
mogeneity scales / = kL with L the background curvature or horizon length scale, the 
following metric provides a consistent description: 

ds 2 = a 2 (r/) [-(1 + 2<p)di] 2 + (1 - 2(f) lij dx i dx j ] , (1) 

(Futamase 1989; JLW1). Here a is the FRW expansion factor (to the order required), r\ the 
conformal time, and 7^ the spatial part of the conformally stationary Robertson- Walker 
metric. The conditions are expected to hold cosmologically everywhere far outside the 
radii of black holes and neutron stars. We call the full time dependence and presence of </> 
in the spatial part of the metric the post-Newtonian corrections. 

The Einstein field equations produce the following relation between the scalar har- 
monic modes of 4> and the density contrast A = dp/p, 

3-0(77, q) + (q 2 + 8iva 2 p - 6k) 0(ry, q) = -47ra 2 pA(r/, q), (2) 
a 

(JLW1; JLW2). Here k is the trichotomic FRW curvature parameter, p the unperturbed 
density, and q is the mode variable conjugate to position. No restrictions are made on the 
size of A, i.e. the density field could be nonlinear. 
Solution of this equation gives (JLW2) 



(3) 



(f)(rj,x) = -(47r/3) ( du(a 3 p/d) ! d 3 y G(u,rj,x,y) A(u,y) 

+ J d 3 yG(7] ,r],x,y)(j)(r] ,y), 

with the Green function in the k = case (used here throughout) 

G(u, V ,x,y) = [a(u)/a( V )} [4nC(u, r/)]- 3 / 2 exp{-|f - y\ 2 /4C(u, V )} 

f v (4) 
C(u,7]) = (1/3) / dw{a/a) . 

J u 
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Transformation of the time integration variable to E = \x — y\/2C 1 / 2 (u, rf) reveals 

/poo 
^ya-^rilx-yl- 1 / dE [a 3 pA(£, y)] e~ E \ (5) 
J Eq 

where Eq = \x — y\/2C 1 ' 2 (r]Q, rf). 

In the Newtonian limit of late times (rj ^> 770) and subhorizon scales (k <C 1), Eq — > 
and it is found that 



(P(r],x) = - / d^ya* — — - 



(6) 



(JLW2), the usual result. [We have neglected the second, initial conditions term from (3) 
due to its exponential die off far from the initial hypersurface, but see §3.3.] 

Going beyond JLW2 we adopt a dust background (a = ao?? 2 = 2H ~ 1 r] 2 where H 
is the Hubble constant) and inhomogeneity behavior A (it, y) = [a(u) / ao]~ 3+p to 

yield 



(j)(rj,x) = 



aoV 71 " 



a 



3* f(y) 



d*y 



\x-y\ Je, 



dE 



1 - 3 



f - y\ 

rf 



-E~ 



-E z 



(7) 



In the Newtonian limit this reduces to <p{r), x) = —a ~ 1 [a(r])/ao\ p ~ 1 J d 3 y f(y)/\x—y\. Some 
particular cases of interest arep = (density unevolving in physical coordinates) and p = 1 
(as for the growth of linear density fluctuations A 1): 



(f)(r],x)p = o 



a(rj) 



d 3 y/^ 
\x - y\ 



erfc(£ ) 



0(77 



,x) p=1 = -- [ d s y-^M- (erfc(£ 
«o J \x - y\ { 



1 + 6 



I —* ->|2 

\x — y\ 
r] 2 



(8a) 



with erfc the complementary error function. 

Since the geodesic equations determining light propagation involve derivatives of the 
potential we calculate in comoving coordinates V0 and ViV 



f 00 



a 



p-i 



\x - y\ A 

rj\2 r°° 



1= I dEU p e~ E2 + 6p |f f dEU p - x E- 2 e- E2 + U p E e 

J Eq ^1 J Eq 



(9) 



Eq 

U = 1-3 



\*-y\ 2 E -2 



with special cases 



V(/> p= o 



a{rf) 



d 3 yf(y)^4s 



V</> p=1 = - / d 3 yf(y) 



a 



\x - y\ 

x-y 
\ x — y\ 



erfc(£ ) + —Eqc'^ 

V7T 



1-6 



I —* -»|2 

\x — y\ 



) erfc(£ ) + -j=E^- E ^ 



(10a) 
(106) 
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The tidal field is V^Vj^, necessary for calculations involving shear of a light ray bundle 
and the resulting image distortions. Because of its length we do not show the expression 
for it here but it is interesting to consider the Laplacian V 2 </>. For p = 



V 2 </> = Ana~ l | f(x) - 7T" 3 / 2 J d 3 yjJ^E$e- E2 ° } 



(11) 



where Co = C{t]q, rj). 

The second term thus illustrates the (not linearized) general relativistic correction 
to the Poisson equation (in the fully specified longitudinal gauge), involving a gaussian 
weighting over the extent of the inhomogeneity. In the "Newtonian" limit the fluctuation 
is restricted to regions much smaller than the dispersion (2Co) 1 ^ 2 so its contrast is averaged 
to zero, leaving the Poisson result, while in the opposite ( "superhorizon" fluctuation) limit 
the gaussian becomes a delta function, giving the Laplace equation appropriate for a 
uniform density field. 

3. Sachs- Wolfe effect 

Given the metric (1) and the results (7), (9), (11) one can investigate light propagation 
behavior by writing down the geodesic equation, geodesic deviation, and the beam or 
Raychaudhuri equation. Properties of the photon bundle, such as convergence and shear, 
and applications of the geodesic deviation equation to astrophysical problems such as 
correlations of observables in terms of the density power spectrum are in ongoing research. 
Here we concentrate on the geodesic equation for individual photon four momentum, in 
particular the Sachs- Wolfe effect on the redshift. 

In longitudinal gauge (<7o« = and proportional to 7^; see Bardeen 1980 and 
Kodama & Sasaki 1984 for gauges and gauge invariance) the expression for the frequency 
of a photon emitted at rj e becomes 

k ( V ) = k ( Ve )[a( Ve )/a( V )](l-2 f dufi-vA . (12) 

(Note the more familiar time derivative of 4> occurs in the unfully specified synchronous 
gauge expression, although one could convert the gradient into a time derivative and a 
surface term.) The ratio of expansion factors is the background cosmo logical redshift and 
h is the photon propagation unit vector. The inhomogeneity induced redshift, i.e. the 
Sachs-Wolfe effect, is then 

• i 

d-qn-Vcj) . (13) 



Upon adopting a density field A(u,y) one can compute the gravitational potential 
gradient by (9) and hence obtain the redshift, or equivalently temperature anisotropy in 
the CMB. This, being an observable, is gauge independent. We consider compact density 
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distributions, both static and time dependent, and then a field of inhomogeneities, along 
with the Fourier representation. 



3.1 Static point mass 

To derive analytic estimates of the dependence of the redshift on density parameters we 
begin by adopting the model of a simple point inhomogeneity of mass m at comoving 
position y = (ycos6>, ysin6>, 0). The x-axis is aligned with the light ray (perturbations to 
the photon path produce effects of smaller order than we consider) and y = 1 — (l+z m ) -1 / 2 
where z m is the cosmological redshift of the mass. The density contrast is A (it, a?) = 
mp~ 1 a~ s (u) 5 ( - 3 ' ) (x — y), corresponding to p = 0. Such a static case is intended as a toy 
model only and is not realistic over long times or large scales. 
For analytic computation rewrite (13) with (10a) as 

z = 2- ds (g2 + 62)3/2 (l - Vl - s)~ 2 [erfc(S) + (2/^)Ee~ E ] , (14) 

where the photon is emitted at x e = 1 — r] e = 1 — (1 + z e ) -1 / 2 , y\ = ycos9, s = x — yi, 
and s = corresponds to closest approach at impact parameter b = ysin.9. We have 
suppressed the subscript on E. First consider the Newtonian case, where the bracketed 
quantity is unity. Around s = the integrand is nearly odd so the dominant term from 
the symmetric part of the interval requires expanding (1 — yi — s)~ 2 to obtain an overall 
term even in s. One readily sees that this gives a logarithmic integral ~ lnb -1 . (Although 
the mass is static, the background is not so an energy shift is expected; but see the next 
subsection for the p = 1 case.) 

Another interval giving a potentially large contribution is when s ~ x e — y\ or x ~ x e . 
Here the integral varies as (1 — x e )~ l = (1 + z^ 1 ! 2 . Physically this corresponds to the 
early epoch when the universe was smaller and so the mass was nearer the source, or 
equivalently the mass per physical volume was greater. The final order of magnitude 
estimate is therefore z ~ (m/ao)[(l + Ze) 1 ! 2 + ln& -1 ]. For small angles b ~ yd. In terms 
of our original parametrization m/ao ~ e 2 K and 9 > k so the "Newtonian" result is 
z N < e 2 K[(l + z e y/ 2 + hiK" 1 ]. Recall that m/a w 10- 8 (m/10 15 M o ), 10" w 5 x 10" 5 
radians so \iik~ 1 pa 10, and the microwave background, for example, originates from 
z e « 10 3 so z N = O(100mH < 10" 6 ). 

In the post-Newtonian case, near the mass the integrand is approximately the same, 
the correction factor being 1 + 0(E 3 ) 1 + 0(9 3 ), while far away the Green function cuts 
off the contribution exponentially so zpn ~ mH ln^ -1 only. The top half of Table 1 gives 
z in units of tuHq for the static case, showing that the analytic dependences hold well. 

Table 1: Point Mass Sachs- Wolfe Effect 

(z/mH ) N/PN 

9 = 1" 9 = 10" 9 = 100" 



z m = 1 218/134 192/108 166/82 

z m = 3 527/348 453/275 380/201 

A ~ a 



5 



1.93/-1.51 
-0.14/-3.64 



3.2 Time varying point mass 

We now examine m = mo(a/ao) p , corresponding to an accreting (or evanescing) mass. In 
density perturbation theory inhomogeneities grow due to gravitational instability as a 1 in 
a k = dust background. Here of course we are not restricted to linear fluctuations A <C 1; 
the parametrization, as can be seen from (2) or as discussed in JLW2, is O(A) = e 2 /^ 2 , 
which can be large. 

The analog of (14) in the Newtonian limit is 

/x e —J/i 
^ , 2 ' 3/2 (l-yi-^) 2(p - 1) . (15) 

Again splitting up the interval we find near s = that zn ~ (1 — p)(m/ao)lnb~ 1 . As 
before, the other possibly significant contribution arises near x = x e (for p < 1/2) and is 
found to vary as (l-x e ) 2p_1 = (l+^ e ) (1 ~ 2p)/2 . Thus the total is z N « mH [(l-p) ln6» _1 + 
(1 + z e )( 1-2p )/ 2 ]. This of course agrees when p = with the static result, since 9 ~ b ~ k. 
Inclusion of the post- Newtonian correction again cuts out the 1 + z e contribution. 

Note that for p > 1/2 the logarithm dominates. Also note that for p = 1, as in 
linear theory, we find zn = 2(m/ao)[(s 2 + 6 2 ) -1 / 2 — (s 2 + 6 2 ) -1 / 2 ] so if s a = —s e , i.e. the 
observer and emitter are situated symmetrically with respect to the inhomogeneity and so 
at the same value of gravitational potential, then there is zero redshift. Classically this 
corresponds to a ball gaining kinetic energy as it rolls into a dip then losing it coming 
out. If it rises to the same height on the far side and if the dip is time independent then 
its final velocity equals its initial. We see that for p = 1 the Newtonian potential is time 
independent and only the endpoints contribute - the usual result. 

However, even for p = 1 and symmetric observer-emitter geometry there is a redshift 
(actually a blueshift) in the post- Newtonian case due to symmetry breaking by the distance 
dependence of the correction, i.e. roughly <p ~ a°erfc(i?o) which is neither symmetric nor 
time independent. The shift is of order tuHq, just as for the asymmetric endpoint case. 
The Sachs- Wolfe effect for point masses varying with time as a is given in the bottom 
half of Table 1. Because of the vanishing of the logarithmic term the results are angle 
independent (for small angles) . For p > 1 the logarithmic term reappears and the results 
regain the angle dependence of the p = case. The generalization to an extended matter 
distribution, e.g. a spherical density inhomogeneity, does not significantly affect the results 
(Kendall 1993). 

The time dependent effect for a dynamically evolving isolated Newtonian density 
fluctuation, e.g. a cluster decoupled from the universal expansion, is sometimes known as 
the Rees-Sciama effect (Rees and Sciama 1968). A gravitational dynamical time scale is 
(pA) -1 / 2 ~ n/e so we expect z ~ J di]d t (p ~ ^(pA) 1 / 2 ^, parametrized as k ■ e/n ■ e 2 ~ e 3 . 
Alternately consider an effective mass within a fixed comoving radius at the times a light 
ray enters and leaves the potential well: m(t out ) ~ m(ti n ) [1 + 0(v/c)] with v the typical 
matter velocity, due to infall or peculiar motions. With an impulse approximation (15) 
becomes z ~ mHQb~ l v ~ e 3 since v ~ (p 1 / 2 ~ e for bound systems. From this order 
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of magnitude parametrization it is clear that to obtain a frequency shift (or equivalently 
temperature shift in the CMB) large enough to be observable (say > 10 -6 ) one needs 
a rapid variation in the potential, more rapid than the gravitational time scale, such as 
from relativistic cosmic strings or black holes (but recall that the metric (1) was derived 
for nonrelativistic, weak field perturbations). Thus, within the context of a dynamical 
Sachs- Wolfe effect, the post- Newtonian formalism provides a significant but unobservable 
difference. 

3.3 Fourier decomposition 

For a more diffuse density distribution, e.g. a linear density field, it is convenient to 
Fourier decompose the gravitational potential into modes corresponding to a characteristic 
fluctuation wavelength, or density inhomogeneity length or mass scale. In fact, working 
in Fourier space allows simplification of many expressions. For example, to obtain the 
Laplacian V 2 </> for arbitrary p, just transform (3), multiply by g 2 , and perform the inverse 
transform. One finds that for general p, (11) is altered by a term (a/ao) p multiplying 
f(x) and a sum inside the integral of polynomials up to order p involving \x — y\ 2 and Co, 
multiplying the exponential. 

Having seen that the variation of the gravitational potential over the photon propa- 
gation gives an insufficiently large effect to be readily observable, we are led to consider 
the effect of the potential at the endpoints of the path, i.e. at emitter and observer, on the 
photon energy. This gravitational redshift is simply z = e — O . Whereas in the isolated 
mass cases of §3.1, 3.2 we could ignore the initial condition term in (3), either by saying 
it provided a term in the redshift that just added independently to the propagation effect, 
or by realizing that its contribution was exponentially suppressed far from the mass, now 
when dealing with a density field and an endpoint effect we must treat it more carefully. 
Either we can ignore it by pushing i] — > 0, or include it in the calculation. We choose the 
latter. 

We concentrate on the last scattering surface of the microwave background radiation 
and the density perturbations there, writing the wavenumber from now on as k instead 
of q to agree with the usual notation in the literature. In the linear density perturbation 
regime A grows linearly with scale factor a, corresponding to our case p = 1. The mean 
square gravitational redshift is proportional to the sum over Fourier modes of the square of 
the gravitational potential, assuming random phases between modes: (z 2 ) ~ J d 3 k \4>k\ 2 - 
The magnitude of the temperature anisotropies (AT/T) 2 = (z 2 ) is usually characterized by 
the zero separation correlation function C(0) but our case is slightly different. For lines of 
sight separated by some large angle ip the two point correlation function C(ip) = (z(0) z{i^)) 
does not vanish because of the coherence introduced by the (4> 2 ) term. Rather than have 
(AT/T) 2 = 2[C(0) - C(V0] -> 2C(0) for V large, one has (AT/T) 2 -> (</> 2 ), i.e. (<j) 2 e ) plays 
the role of the zero lag correlation in setting the magnitude. Alternately one can say that 
<po offers only a constant, isotropic shift so one can ignore it and consider only z = <p e . 

This quantity can be written in terms of the density perturbation power spectrum by 
taking the Fourier transform of (8b), but first we must evaluate the initial condition term 
of (3) and add it in. It is 

( p i : c) (r ] ) = la(r ] o)/a(v)]Mrio)e- k2c °. (16) 
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Note the exponential suppression as r\ grows larger (later) than 770- At this point we see 
that we really have two inputs to specify in (3), the initial density field A(??o) and the 
initial potential (p(r]o). Because of the form of the full relativistic equation (2), both must 
be given. This is not unexpected though, because in JLW2 it was pointed out that (2) 
was essentially analogous to an inhomogeneous diffusion equation with time dependent 
parameters. It was found that the potential corresponded to the temperature in that 
sort of problem, and the density contrast to a heat source term. From that physical 
situation, however, we know that we must generally specify the initial distribution of both 
the temperature and the heat sources. It is only in the late time limit that the source 
completely determines the diffusive variable; in our case this corresponds (if compactness 
holds as well) to the Newtonian limit and the Poisson-Newton equation. 

For the initial condition constant 4>k(vo) we use the ansatz 4>k(vo) = —^Ak~ 2 fk, so 
A gives the deviation from the Poisson equation (A = 1) in the initial conditions. Adding 
(16) to the transform of (8b) gives in terms of the density perturbation power spectrum 

(Z 2 ) = 167T 2 f d 3 kP k k- 4 C k , 

C fc = [l - (1 - A)( m /v) 2 e- k2Co - 12fc" V 2 (l " e" fc2Co )] , 

where Pk = |Afc| 2 is the density power spectrum and Cq — C(rjo, rj) as given by (4). 

Note that C gives the overall correction factor to the usual Poisson relation (in Fourier 
space) between the gravitational potential cf> entering the metric (1) and the energy density; 
in particular, C = A 2 at rj = tjq. It can also be obtained directly from transforming the 
p = 1 version of the Laplacian (11). Thus C can be viewed alternately as giving the post- 
Newtonian adaptation of the endpoint Sachs- Wolfe effect, or as an effective alteration to 
the density power spectrum Pk, i.e. an effective transfer function. It is included in any 
(e.g. numerical) treatment using the full equation (2). 

Depending on the initial conditions it has the possibility of either enhancing or sup- 
pressing the low Fourier modes, i.e. increasing or decreasing the large scale power in the 
intrinsic power spectrum This, for example, would cause a smaller (larger) overall 
normalization factor to be needed to match the large angle COBE microwave background 
anisotropy measurements, and hence also decrease (increase) the resultant predicted small 
scale power, thus ameliorating (exacerbating) the difficulties of the cold dark matter model. 

The correction factor Ck is plotted in Figure 1. In the late time limit, rj ^> 770, the 
results are independent of 77/770 and A. From the expression in (17) or the Figure, three 

— 1/2 

regimes in perturbation wavelength can be identified. When k ^> C , the factor is 
close to one, i.e. this is the late time, compact inhomogeneity Newtonian limit. For 
\[Ylr)~ x Cq 1 ^ 2 , the factor begins to decline predominately due to the influence 

of the exponential, and when k <C \J~Ylr}~ 1 , the {krf)~ 2 term becomes important as well, 
leading to Ck — > A 2 (rjo / rj) 4 . From Figure 1, the deviation from the usual Pkk~ 4 behavior 
of the Sachs- Wolfe effect becomes noticeable below kr) = 10 — 40, depending on the value 
of Co- Converting to wavelengths by 

kr] = 47tri/(H \o) = 120 (Xo/lOh^Mpc)' 1 [(1 + z e )/10 3 ]" 1/2 , (18) 
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with Ao the present day wavelength and z e the redshift of the last scattering surface, we 
see this corresponds to scales A > 30 — 100 h~ x Mpc or angular scales 9 > (1/2) — 2°, 
applicable to the COBE regime. 

4. Conclusion 

The pseudo-Newtonian gravitational potential appearing in the metric (1) that realistically 
approximates our universe, incorporating inhomogeneities in a FRW background, is related 
nontrivially to that density fluctuation distribution. The Green function (or kernel or 
propagator, in field theory terms) generally has a different spatial extent, or compactness, 
than the Newtonian case which corresponds to its late time, near neighborhood limit. In 
consequence, light propagation behavior for those regions of parameter space far from that 
limit differ from the Newtonian. 

In the dynamical Sachs- Wolfe effect of the redshift induced by inhomogeneities along 
the light path, the deviation was shown to be significant but the overall magnitude too 
small to be observationally interesting. For the endpoint Sachs- Wolfe effect of gravitational 
potential fluctuations at the source, the deviation could be viewed as an effective transfer 
function modifying the density perturbation power spectrum, with potentially important 
consequences for the normalization of the primordial density spectrum and hence the 
predictions of that model on smaller scales as a galaxy formation scenario. 

Acknowledgments. I thank Mark Jacobs, Jason Kendall, and Bob Wagoner for useful 
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Figure caption 

Fig. 1. The Poisson correction factor Ck = |</>fc| 2 /(167r 2 /c~ Pk) is plotted vs. krj for various 
initial conditions (A, 77/770). For small wavenumbers (large scales) a significant enhance- 
ment or suppression can be achieved. The late time curve 77 ^> 770 is independent of A. 
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